Monovalent Ion Condensation at the Electrified Liquid/Liquid Interface 
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X-ray reflectivity studies demonstrate the condensation of a monovalent ion at the electrified interface between 
electrolyte solutions of water and 1,2-dichloroethane. Predictions of the ion distributions by standard Poisson- 
Boltzmann (Gouy- Chapman) theory are inconsistent with these data at higher applied interfacial electric 
potentials. Calculations from a Poisson-Boltzmann equation that incorporates a non-monotonic ion-specific 
potential of mean force are in good agreement with the data. 
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Interfacial ion distributions underlie numerous electro- 
chemical and biological processes, including electron and 
ion transfer across charged biomembranes and energy 
storage in electrochemical capacitors. The solution to 
the Poisson-Boltzmann equation for a planar geometry, 
Gouy- Chapman theory, including modifications with a 
Stern layer, is often used to predict ion distributions 
near those interfaces.^ We showed previously that 
the predictions of such theories are inconsistent with 
x-ray reflectivity measurements of ion distributions 
at an electrified liquid/liquid interface.^ Instead, an 
ion-specific Poisson-Boltzmann equation (PB-PMF) 
that incorporated a potential of mean force (PMF) for 
each ion produced excellent agreement with the x-ray 
results. 1^ The PB-PMF theory accounts for interactions 
and correlations between ions and solvents that are 
left out of Gouy- Chapman theory. This approach has 
promise for understanding ion-specific effects that are 
relevant to many chemical processes.^ 

Here, we demonstrate the condensation of a mono- 
valent ion at a liquid/liquid interface. Recent theory 
proposes that condensation of multivalent ions is the 
result of strong ion-ion correlations.^ However, these 
theories do not predict such distributions for monovalent 
ions. Our current results can be understood by PB-PMF 
theory. We have chosen to fit our x-ray data to the 
potentials of mean force instead of fitting to a model 
of the electron density profile because a single PMF 
for each ion determines the ion distributions for all 
interfacial potentials. 

The system under study is the liquid/liquid interface 
between a 100 mM aqueous solution of NaCl (Fisher 
Scientific) including 20 mM HEPES to buffer the pH 



to 7.0, and a 5 mM solution of bis(triphenyl phos- 
phoranylidene) ammonium tetrakis(pentafluorophenyl) 
borate (BTPPA+, TPFB") in 1,2-dichloroethane (DCE, 
Fluka). Water was produced by a Barnstead Nanop- 
ure system and DCE was purified using a column of 
basic alumina. BTPPATPFB was synthesized from 
BTPPACl (Aldrich) and LiTPFB (Boulder Scientific). ^ 
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FIG. 1. Circular glass sample cell and x-ray kinematics. Elec- 
trochemical cell diagram: Ag|AgCl | O.IM NaCl | water | +20 
mM HEPES | 5 mM BTPPATPFB | DCE | 10 mM LiCl+1 
mM BTPPACl | water | AgCl | Ag. A four-electrode poten- 
tiostat (Solartron 1287) is used to apply potential at counter 
electrodes (CEi,2, 9 cm^ Pt mesh) and monitor potential at 
reference electrodes (REi,2) in Luggin capillaries within 4 mm 
of interface. The liquid/ liquid interface of 7 cm diameter 
is pinned by a Teflon strip (affixed to the glass wall) and 
flattened by adjusting the volume of DCE phase. Volume 
ratio of water: DCE is 2:1. The x-ray wave vector transfer 
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Conductance measurements using the method in RefP 
determined that 54% of BTPPATPFB is dissociated in 
DCE. 

The electric potential difference A(/)^~^(= 
^water _ ^oii^^ between the water and oil (DCE) 
phases is given by the applied potential difference across 
the electrochemical cell (FigfTl) minus the potential 
of zero charge ° = A^^^",, - A0»-°). We 

determined A0^^~^ = 318 ± 3 mV by measuring the 
interfacial tension as a function of ^^i^appUed^' ^^^^ 
Na+ and Cl~ stay primarily in the aqueous phase and 
BTPPA+ and TPFR- stay in the DCE phase through- 
out the potential range studied. When A<p'^~^ ^ 0, 
the ions form back-to-back electrical double layers at 
the interface. For example, when A^^~^ > 0, the 
concentrations of Na+and TPFB~ are enhanced at the 
interface while those of Cl~ and BTPPA+ are depleted. 
The variation of ionic concentration along the interfacial 
normal produces a variation in the electron density 
profile p{z) (averaged over the x-y plane) that is probed 
by x-ray reflectivity. 
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FIG. 2. X-ray reflectivity R{Qz) normalized to Fresnel re- 
flectivity Rf{Qz) for various potentials across the water/ 1,2- 
dichloroethane interface as a function of wave vector transfer 
normal to the interface at T=296 K. From top to bottom 
A0"^-" = 0.33 V (o), 0.28 V (•), 0.18 V(o), 0.08 V (•), -0.02 
V(o), and —0.12 V(o) (offset for viewing purposes). Dashed 
lines: Gouy-Chapman theory. Solid lines: PB-PMF. Red and 
blue lines indicate the use of two different PMFs for TPFB~ 
(see text). 



X-ray reflectivity measurements R{Qz) from the 
electrified liquid/liquid interface were carried out at 
the ChemMatCARS sector of the Advanced Photon 
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FIG. 3. Potentials of mean force for BTPPA+ (black) and 
TPFB" [W^pp^-{z):Yed,W^pp^-(z)-h\ue] determined by 
fitting the reflectivity data in Fig. l2| PMFs for Na^ (green 
dots) and Cl~ (circles) were calculated by MD simulations 
(see text) (Ref.™. 



Source.^ R{Qz) is the reflected intensity normalized by 
the incident intensity (after subtraction of background 
scattering^) as a function of wave vector transfer Qz = 
(47r/A)sina , where A (=0.41255 ±0.00005 A) is the x- 
ray wavelength and a is the angle of incidence (Fig. 1). 
Figure 2 illustrates R{Qz)/Rf{Qz) for different A0^-^. 
The variation of the peak amplitude in R/Rp with in- 
creasing A^^~^ reveals the formation of a TPFB~ layer, 
as discussed below. 

We analyzed the x-ray data using the Poisson- 
Boltzmann (PB) equation, 

^ i 

which relates the electric potential (p{z) along the in- 
terfacial normal z to the concentration profile of ion i, 
Ci{z) = ex.-p[— AEi{z) / ksT] , with Boltzmann constant 
ks, temperature T, charge of ion i (BTPPA+, TPFB", 
Na+, and Cl~), permittivity of free space ^o, and dielec- 
tric constant s of either DCE (10.43) for z < or water 
(78.54) for z > . AEi{z) is the energy of ion i relative 
to its value in the bulk phase. The bulk ion concentration 

is calculated from the Nernst equatioiff^ Fitting 
to the data involves calculating the electron density p{z) 
and R/Rf from the ion concentration profiles Ci{z) as de- 
scribed previously.^ For the purpose of calculating p{z) 
from Ci{z) the ions were modeled as spheres of diameter 
2 A for Na+, 3.5 A for Cr, 12.6A for BTPPA+, and lOA 
for TPFB~, where the latter were estimated from the 
crystal structure of BTPPATPFB.^^ The TPFB" ion 
provides the dominant ionic contribution to the electron 
density profile when A^^~^ > . 

The Gouy-Chapman theory assumes that Ei{z) = 
ei(j){z) in Eq. [l] Fits of R/Rp to predictions of Gouy- 
Chapman theory (Fig. |2j dashed lines) used only the in- 
terfacial roughness and a Qz offset ( 10~^A~^, a typical 
misalignment of the reflectometer) as fitting parameters. 
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TABLE I. TPFB" PMF parameters obtained by fitting to the reflectivity data at A^™"" = 0.28 and 0.33 V. 
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-5 ±0.5 
-25 ±0.5 


3 ±0.1 
11 ±0.5 


9±4 
10 ±2.7 


-3.5 ±0.2 
-7.5 ±0.3 


3.4 ±0.2 
2.6 ±0.3 


-9 ±0.25 
-5.25 ±0.2 



These fits agree with the data at small A^^ ^ (—0.12 
V to 0.18 V), but at larger A^^-^ (0.28 V and 0.33 V) 
R/Rf is greatly overestimated primarily because Gouy- 
Chapman theory predicts unphysically large ion concen- 
trations near the interface. 

Ion-specific effects can be included in Eq. [l]by express- 
ing Ei{z) ~ ei(j){z) ± Wi{z)^ where Wi(z) is the potential 
of mean force (PMF) for each ion i.EEMU xhe PMF of 
Na+ was calculated from a molecular dynamics (MD) 
simulation for a single ion.^^ The PMF of Cl~ was taken 
from an MD simulation in the literature.'^ Fig. [s] illus- 
trates the monotonic variation of Wi{z) for Na+ and Cl~. 
Due to the computational difficulties of simulating Wi{z) 
for large molecular ions such as BTPPA+ and TPFB~ we 
used a phenomenological PMF previously introduced in 
ReP, 

where p{= w^o) refers to either the water phase {z > 0) 
or the oil phase (DCE,z < 0), - is the Gibbs en- 
ergy of transfer of ion i from water to oil, is an offset to 
ensure continuity of Wi{z) at z = 0, and characterizes 
the decay of Wi{z = 0) to its bulk values W^"" andVI^f . We 
used this monotonic PMF for BTPPA+, but had to mod- 
ify it for TPFB-, as described below. Since - 
for BTPPA+ is known (Ref. ,the PMF of BTPPA+ 
is characterized by 3 parameters determined by fitting to 
R/Rf data at A^^"^ = -0.12 V where it is expected 
that the BTPPA+ interfacial concentration is enhanced: 
^BTPPA+(= 14±12/-6A), L^j,FPPA+i= 20±11/-6A), 
and Wbtppa+{^){= 13 ± 2A:bT). The large error bars 
on the PMF of BTPPA+ are due to the small magnitude 
of the most negative A(/)^~^ that we studied. 

The x-ray refiectivity at the two highest positive poten- 
tials cannot be fit if Eq. [2] is used to model the PMF for 
TPFB~. The simplest model that will produce the peaks 
in Fig. [3] is a single layer of TPFB~ ions at the interface 
(note that a layer of Na+, whose concentration is also en- 
hanced at the interface, cannot provide the x-ray contrast 
required to fit the data). The TPFB~ layer is modeled 
by an attractive well in the PMF. W^pfb- (^) is given by 
Eq. [2]plus a Gaussian function D ex.p[—{z — zo)'^ /2apj^j^] 
for z < along with a constant offset at z = to main- 
tain continuity (see Figj3|. Analysis with a Lorentzian 
produced similar results. 

The six parameters of Wtpfb- {^) [^o, ctpmf^ 



^TPFB- ' ^TPFB- ' ^TPFB- (0)] along with the Qz offset 
and the interfacial roughness ( 4. 3 A < a < 5.1 A) are de- 
termined by fitting R/Rp measured at A0^~^ = 0.28 
V and 0.33 V, where the concentration of TPFB~ is 
enhanced at the interface (Table |l| . This fitting is per- 
formed under the constraint that the resultant Wion(^) 
produces R/Rp in agreement with the data over the 
entire range of potentials. In addition, fitted PMFs 
were rejected if the fit value of the roughness <j was un- 
physically small. In those cases an interfacial bending 
modulus^^ on the order of IOOO/cbT would have been re- 
quired to reconcile the discrepancy of a with its value 
predicted by capillary wave theorj/^^. In the case of the 
TPFB~ PMF, two local minima in x^-space (denoted 
^TPFB-(^) ^TPFP-(^) ) were found to satisfy 
these conditions. Potential profiles that are intermediate 

between W^p^^^ {z) and w![Jppp_ do not satisfy these 
conditions. Most of these fits had values of a within 
one standard deviation of capillary wave theory predic- 
tions using the measured potential-dependent interfacial 
tension."^. Fits to W^^ppp_ {z) at A^^-^ = 0.28 V and 
0.33 V had values of a within two standard deviations of 
capillary wave theory. 

The PB-PMF model with the Wi{z) shown in Fig. [s] 
produces R/Rp in good agreement with the data over the 
entire range of measured potentials (Fig. |2]). The attrac- 
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FIG. 4. (left) Ion concentration profiles (in units of molar- 
ity) at A0^-^ = 0.33V calculated from PB-PMF. BTPPA+ 
(black), TPFB- [W^pFB-(^) :red,M^^^^^_(z) :blue], Na+ 
(dots) and CP (circles), (right) Electron density profiles for 
various potentials calculated from PB-PMF. Top to bottom: 
A(/)"^-" = 0.33V [W^pp^_{z):Ted,W^^pp^_{z):Uue], 0.28V 
(dashed), 0.18V (solid), 0.08(dashed), -0.02V (solid), and 
-0.12 V (dashed). 
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tive wells for Wj.pj^^_ have comparable depths {6 ksT 
for W^pj,^_{z) and 5 ksT for V^^^^^_(z)), FWHM, 
and centers (Table 1). The ion concentration profiles 
are calculated from Eq. [l] using Wi{z). Figure |4] 
shows that the Ci{z) at the highest potential, A^^~^ = 
0.33 V, take the form of two back-to-back double layers 
with a sharply defined layer of TPFB~. The different 
Ci{z) calculated from or W^^^^_ (z) differ 

mainly in the broadness of the profile, which in the case 
of W^pjpp_ (z) returns to its bulk value at z = 0, while 
^TPFB- (^) allows TPFB~ to penetrate slightly more 
into the water phase. The electron density profiles p{z) 
calculated from the different Ci{z) are almost identical 
(Fig. |3|, which demonstrates why our data cannot dis- 
criminate between Wj.ppp_ (z) and W^ppp_ (z) . 

The maximum density of TPFB~ near the interface 
occurs at A^^-^ = 0.33 V and is 1 nm^ per TPFB" 
ion when Wrj.ppp_ (z) is used or 1.5 nm^ per TPFB~ ion 
when W^ppp_ (z) is used. Both values represent a high- 
density layer for an ion of 1 nm diameter. Although dense 
ionic layers have been observed in the interfacial adsorp- 
tion of charged amphiphiles,^^ the absence of a dense 
TPFB- layer at A0^-^ ^ indicates that TPFB" is, at 
most, weakly amphiphilic. 

Simulations and supporting spectroscopy experiments 
indicate that highly polarizable ions (such as I~ with a 
polarizability of 7.4A^) are preferentially adsorbed to the 
water/vapor interface, though dense layers are not ex- 
pected or observed. We calculated the polarizability of 
TPFB- to be 42.91^.22^^^ This large polarizability may 
play a role in forming a dense layer at high potentials. 
Also, Borukhov et al suggested that the entropy of the 
solvent can stabilize large ion adsorption. Additional 
theoretical work is required to determine the relevance 
of these two effects for the data presented here. 

The MD simulations of the potentials of mean force 
that we used for Na+ and CI" do not account for ion-ion 
correlations, but they do include ion- solvent and solvent- 
solvent correlations. Such correlations also account for 
the monotonic form of Wbtppa+{^) • However, as a 
result of modeling the x-ray refiectivity, the phenomeno- 
logical Wtpfb- {^) Fig. [s] must implicitly account 
for ion-ion correlations if they are important for the ob- 
served condensation. The description of this monovalent 
ion condensation within PB-PMF theory illustrates the 
utility of this approach in describing ion-specific effects 
that are important for the behavior of ions in soft matter. 
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